Particle motion in weak relativistic gravitational fields 
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We derive the geodesic equation of motion in the presence of weak gravitational fields produced 
by relativistic sources such as cosmic strings, decomposed into scalar, vector and tensor parts. We 
find that the vector (gravito-magnetic) force is an important contributor, and for non-relativistic 
particles we recover the well-known result for the impulse from a moving straight string. Our results 
can be straightforwardly incorporated into N-body simulations to allow for the presence of cosmic 
defects or other sources of weak gravitational fields. 
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I. INTRODUCTION 

Topological defects such as cosmic strings [H-Q 
generic by-products of many inflationary models 
and of Grand Unification adding a characteristic sig- 
nature to the gravitational and matter fiuctuations pre- 
dicted by inflation. Precision Cosmic Microwave Back- 
ground (CMB) measurements set limits on the allowed 
defect abundance and thus also on infiationary 

models that produce defects. 

There is also accurate data on the galaxy power spec- 
trum over a wide range of scales [l6j , which constrain the 
matter perturbations. Inflation creates a nearly Gaus- 
sian spectrum of initial perturbations (see e.g. [l7|). 
whose subsequent evolution can be computed in linear 
theory, and compared to the data (under assumptions 
about the bias, i.e. the ratio between the galaxy and 
matter power spectra). However, defects are localised 
and "active" [l^ sources of gravitational fields, creating 
highly non-Gaussian perturbations [l9l - [2ll |. For example, 
cosmic strings create a wake behind them as they move 
through matter [22l . [23| , in which there is a planar rela- 
tive overdensity of order 1 as soon as it is created. The 
evolution of the wake is therefore immediately non-linear, 
and standard linear theory in Fourier space does not ap- 
ply. The strong non-Gaussianity is very likely to impact 
the growth of structure and (for example) could affect 
the bias of the galaxy power spectrum. For this reason, 
we need to find ways to go beyond linear perturbation 
theory to calculate the matter power spectrum derived 
from the gravitational perturbations of defects. 

One way to do so involves N-body simulations. In 
some early numerical work by (23 . [2a | the structure of 
the wake induced by a single straight string was studied 
by setting up a velocity kick towards the plane behind 
the string as the initial velocity perturbation, as derived 
from semi-analytical predictions |l|, [2^ , verifying the pre- 
dicted width of the wake, the inflow velocity of the dark 
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matter and assumptions about the self-similarity of the 
solution. 



However, in general defects have a complex and evolv- 
ing 3-dimensional structure that extends up to the hori- 
zon scale. For example, strings are not straight but form 
a tangled, self-intersecting and fast-moving network of in- 
finite and closed pieces, with a characteristic length scale 
of about one third of the horizon scale and a character- 
istic speed of about a half that of light [13, [ll] . 



In order to capture the full impact of the defect per- 
turbations on the large-scale structure, we will need to 
work on cosmological scales and with relativistic sources. 
For this we will need general relativity, both for the per- 
turbations in the gravitational field, and the deviations 
to the motion that the field produces. 



In this paper we derive from first principles the 
equation of motion of massive particles in a perturbed 
Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmol- 
ogy, keeping all terms linear in the gravitational fields, in- 
cluding vector and tensor. This is required for relativistic 
sources such as topological defects, where all parts of the 
energy-momentum tensor are comparable in magnitude. 
We will then test our formalism with the help of a moving 
straight string, for which the impulse on passing particles 
is known exactly. Finally, we perform an N-body simula- 
tion to study the growth of the wake behind the moving 
string, comparing with previous work using Newtonian 
gravity only [2^ . 



This paper establishes the formalism by which the ef- 
fects of topological defects on the matter in the universe 
can be taken into account, as a necessary preliminary to 
non-perturbative (especially N-body) calculations of the 
growth of structure in cosmological models with defects. 
The formalism is more general, however, and allows to 
add any sources of weak relativistic gravitational fields 
to N-body simulations. 
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II. OVERVIEW OF LINEAR PERTURBATION 
THEORY 

Very generally, given a metric g^i, we can compute the 
Christoffel symbols 

^ap = ^9'"' {dagpy + dpgau - d^g^p) , (1) 



the curvature tensor i?^^ and thus the Einstein tensor 
G^i, as well as the geodesic equation of motion, 



(fx'' 



0. 



The gravitational field equations or Einstein equations. 



(3) 



describe the interaction between gravity and matter, with 
the latter given by its energy- momentum tensor T^^,. 

In general the Einstein equations are very difficult to 
solve, but in the limit of weak fields and small pertur- 
bations we can use linear perturbation theory around a 
fixed background metric |29l - [3]| . These linearised equa- 
tions then naturally decompose into irreducible compo- 
nents under rotations, scalars (S), vectors (V) and tensors 
(T). However, at higher order in perturbation theory this 
decomposition is no longer maintained and the different 
types of perturbations mix. 

We will assume throughout this paper that the sources 
of the gravitational perturbations evolve on the back- 
ground universe and are not affected by the perturba- 
tions that they generate, i.e. that their perturbations 
only affect the remaining constituents. This is the case 
for topological defects, for which we can perform the 
numerical simulations separately, recording the gravita- 
tional perturbations which are then self-consistently in- 
serted into the linearised Einstein equations of the full 



system [32|. By running the defect and the N-body sim- 
ulation in parallel and exchanging information between 
the two we could in principle include the so-called grav- 
itational backreaction on the defects, but this is left for 
later work. 



A. Metric and Christoffel symbols 

We choose a background metric a^rj^i,, where rj^i, = 
diag(— 1, 1, 1, 1) is the Minkowski metric (i.e. we only 
(2) consider flat space) and a perturbation a^hi^^y so that the 



full metric with the S,V,T decomposition becomes 



(4) 



The scalar perturbations in the conformal Newtonian 
gauge, the vector perturbations in the vector gauge and 
the gauge invariant tensor perturbations are defined by 



hlfjdx^dx" 



h^^^Uxt'dx" 



-2iljdT^ + 2(l)5rjdx'dx^ (5) 
-2Y.,dTdx' (6) 



fc'S,- = 



h'-^Jdx^'dx" 



hlPdx'dx^ 



(7) 







where r is conformal time. 

From the metric we can immediately derive the 
Christoffel symbols with the help of Equation ([Ij. Us- 
ing primes (' = dr) to denote derivatives with respect to 
conformal time, we find 
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In particular, the zcroth component is found to be 



AT) 



AT). 



and the i'th component is found to be 



(8) 



(9) 
(10) 
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B. Einstein Equations 

In Fourier space, and following the formalism of 
the Einstein equations can be written as 

a 



(11) 



iP = -SttGU 



16nG 



,iv) 



n 'J 



where we used the following elements of the energy mo- 

(12) 



mcntum tensor T^i,: 




— Too 


fv 




fp 


" 3 






f. 




(V) 


= {Toi - hk^Toj) 




= {PaP,rn-{ll2)P^,Phn)P'^''P'''Tab 



where hats denote unit vectors, and where we used the 
projection operator 



Pij — ^ij kikj . (13) 



III. PARTICLE MOTION IN WEAK 
GRAVITATIONAL FIELDS 



Evaluating the zeroth component of the geodesic equa- 
tion ([2]) gives us an expression for the evolution of the 
energy of massive particles 



(1 - 2'ijj + 24>)xjx^ - 2Y.,£ + hiPx'x^ 



(14) 



Using this in the i equation we find the equation of 
motion of massive particles in a weak gravitational field: 



ip-2<j> + 2dj{'4} - (j))x^ + (j)XjX^ + {djEk + iA|p)f '^f^' 

2 



-1 + (1 - 2i/) + 2(t))xjX^ - 2T.jX^ + h'j^-^xkx^ 



~di{i}} — (pXjX-') 



djTiiX-' — di'SjX-' 



dkh 



^PxH*' + IdJi'iJ^M'^x^ . 



kj 
d 
a 



^i{XjX'^ 



r 



(15) 



We notice that all types of perturbations, scalar, vec- 
tor and tensor, affect the particle motion. However, to 
leading order for non-rclativistic particles with x 1 we 
find that tensor perturbations do not contribute. 

Converting to physical time and retaining only the 
leading order terms in the particle velocity we find 



x' + 2-x' 
a 



1 



1. 



— E 



(16) 



potential of the of the topological defects.^ However, we 
find that vector-type metric perturbations, if they are 
present, affect particle motion at the same level. Thus 
gravito-magnetic forces can be as important as the stan- 
dard scalar force for a relativistic source. This expression 
should be used when extending N-body codes to account 
for general sources of weak gravitational fields. The full 
result needs to be evolved if particles can reach relativis- 
tic speeds, and it may be worth to occasionally monitor 
the size of the next order terms. 



The second term on the left is due to the expansion of 
space. Although it contains x it needs to be taken into 
account as it is not suppressed by any of the metric per- 
turbations. The first term on the right is the usual gradi- 
ent of the gravitational potential which is the sum of the 
potential due to the particles and the scalar perturbation 



^ We allow for a non-zero anisotropic stress since for topological 
defects in general 7^ -0. Wc note that only the ^ potential ac- 
celerates massive particles to lowest order in the particle velocity. 
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IV. PARTICLE MOTION INDUCED BY A 
STRAIGHT STRING IN MINKOWSKI SPACE 




FIG. 1: A Nambu Goto string is aligned with the z axis travel- 
ling in the x direction at constant speed v through the middle 
of the xy plane. A test particle is at a position ( Aa;, Ay) w.r.t. 
the initial position of the string such that Aj//Ax <C 1. The 
string travels a distance 2 As. 

In this section we will reproduce the well known ef- 
fect of a velocity kick induced by a moving straight 
Nambu-Goto string on a test particle, in a non-expanding 
Minkowski background with a = 1. We compute the 
result by evolving equation ([T6| . which illustrates the 
approach that we will be using in the future for the full 
network and verifies that we arrive at the correct answer. 
Wc show that both scalar and vector perturbations con- 
tribute significantly to the particle motion, even though 
in the approximation where the particle is scattered by a 
string coming from infinity, in the infinite time limit only 
the scalar contribution remains relevant. 

The string Nambu Goto action is given by 



S 



(17) 



where /i is the string mass density, (a = 0, 1) are 
the coordinates on the worldsheet traced out by the 
string with spacetime coordinates X^(C), and g'^' is the 
determinant of the induced metric on the world-sheet. 

The particular gauge chosen for defining the world- 
sheet is PI 



(2) 
/Ol 



0, 



(2) 

ffoo 



(18) 



and we may also identify worldsheet and coordinate time 
with the choice = t. 

Wc are making the simplification of the string living 
in an flat, non-expanding space. Its energy momentum 
tensor is [i[ 

T^'''[x, t)=iij dX (x^X^^ - l^l"^) 5^3) _ x{X, t)) 

(19) 

where = ^qX^ and X^' = 



We are considering specifically a straight string parallel 
to the z axis, traveling in the x direction at a constant 
velocity w ^ so that 



(l,i>,0,0), XA' = (0,0,0,i), (20) 



In this case we find that the energy momentum tensor is 
given by 

T'^-'ix, t) = ^fiAP^'^Six -xo- vt)S{y - yo) (21) 



with 



f 1 V 
V 




(22) 



Vo -1/7V 



We now change to Fourier space in order to solve for 
the perturbations in the metric with the help of the Ein- 
stein equations: 

T'"'{k, t) = 2TTjnM^"'e'''^^''°+'"'^e'''yy"S{k^) (23) 

Using equations (|lll) with a = 1 we find the scalar and 
the vector perturbations to be 



(j){k,t) 



,ik^{xo+vt) iky-ya 



(24) 



- 2v' 



327r2GAi7 



fc2 



Vi kihxV 



xe 



ikj,{xo+vt) ikyyo 



SiK). 



(25) 



(26) 



In our setup, we have that Vi = (w, 0,0). To calculate 
the effect of the perturbations on a particle to first or- 
der (equation [T6| we only need tp and Si. Their inverse 
Fourier transform is given by 



2r2 
2Gii'yv 



-Qx^ 



-2x^ 



log -2 



AG^jvxy 







(27) 

(28) 

(29) 
(30) 



where tq is an integration constant that sets the dis- 
tance at which the logarithmic contribution of the in- 
finite straight string to ip and Ei vanishes, and {x,y) = 
{xq + vt,yo)- Redefining the variables 

(31) 
(32) 



X 

y 



Xq + Vt 

yo-y 



^ V in units of c = 1, 7 is the Lorentz factor. 
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equations ([27)) - (|30)) give the fields at position {x,y) 
due to a string at (xg + wt,j/o)- Inserting these expres- 
sions into equation (fT6|) we arrive, after ignoring cohisions 
where the solution is divergent, at the following equation 
of motion for particles: 



ScalcLrField 



Axr,{rl - rl) - 4yry{3rl + r^)] (33) 
A{xry - yr,){rl - rl)] (34) 



We notice that, as expected, the logarithmically diver- 
gent part of the scalar and vector potentials ip a-nd Si 
have been removed by the derivatives and the accelera- 
tion is decaying like 1/r at large distances. The acceler- 
ation field around a straight string is shown in Figure [21 
We solve the equations for x{t) and y{t) numerically 
with a standard numerical solver^ for ordinary differen- 
tial equations and for initial conditions chosen so that 
the string passes close to a test particle in the y direction 
(at a distance Ay), starting at a distance Aa; far away 
and moving past the particle in the x direction again to 
a distance Ax such that Ay/ Ax ^ 1 (the near-by limit). 
We find that the net acceleration in the x direction can- 
cels roughly out, while the net effect in the y direction, 
in particular the velocity kick, agrees with the predicted 
value [H 



(35) 



Our numerical result for the effect of a GUT string with 
G/x « 10~^ travelling at v = 0.333 past a particle with 
Ay /Ax « 0.003 is 



u/up « 0.9985 , 



(36) 



and the agreement can be improved nearly arbitrarily 
as Ay /Ax 0. We also inserted the modified accel- 
eration equations ([33)) and ([34)) into the public N-body 
code Gadget-2 [33] and found that the N-body code re- 
sult agrees with the result from the numerical solver to 
machine accuracy as long as as we resolve the dynami- 
cal time of the particle string interaction in the N-body 
code. This can be tuned to any precision when looking 
at the effect on single particles (however, this may be 
a problem to be solved when doing large scale N-body 
simulations where a general small limit on the maximum 
timestep is too expensive). We conclude that our for- 
malism and simulation set-up reproduces the standard 
results in non-expanding space time accurately. 

Details of the particle motion (from the N-body code) 
are shown in Fig. [3l The acceleration in the x direction 
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FIG. 2: The acceleration field around a straight string in 
Minkowski space-time moving at constant speed in the x di- 
rection (i.e. horizontally through the center of the box). The 
top panel shows only the contribution from the scalar com- 
ponents, and the middle panel from the vector components of 
the metric perturbations induced by the string. The bottom 
panel shows the total acceleration field which is the one seen 
by a particle. 



^ NDSolvc of the Mathematica software paekage. 
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(along the motion of the string, left panels of the figure) 
averages to zero, so that there is no net velocity left after 
the string has passed, and only a small overall displace- 
ment. It is however remarkable how the scalar and vector 
parts combine to a smooth overall motion, which is best 
visible in the middle panel on the left for the velocity 
in the x direction. In the y direction (perpendicular to 
the string motion) we can see the particle receiving the 
above-mentioned velocity kick. The contribution from 
the vector part is small and mostly serves to render the 
kick more step-like. 

When looking in more detail at the late-time impact 
of the scalar and vector parts of the acceleration field, 
we find that in the near-by limit the vector part does not 
contribute significantly to the final velocity, see Fig. 2] 
However, this is due not least to the special case that 
we consider, where a long, straight string moves on a 
straight trajectory past a particle. In reality we will be 
dealing with a string network, in which case strings are 
not straight, and they move on curved trajectories. In 
this case we would not expect to satisfy the near-by limit 
at all times. In this situation the vector part can con- 
tribute at a level comparable to the scalar part. 



When the string passes near a particle, this particle 
will be imparted a velocity towards the string (the veloc- 
ity kick discussed above). In physical coordinates with 
the origin fixed to a point on the string trajectory, how- 
ever, the particle is still following the Hubble fiow until 
the recession velocity drops below the velocity due to the 
string passage. The particle will then start to move to- 
wards the region through which the string has passed, 
and we expect that the particles will form a wake behind 
the string once they reach this region. Quantitatively, 
we can compare the particle motion w.r.t. the axis of 
symmetry and the turn-around radius rt at which the 
particle motion decouples from the Hubble flow to the 
calculation based on the Zel'dovich approximation in [l| : 
The physical particle position y{t) is given by 



y ^a{yo + ^/ai) 



(38) 



where yo is the initial particle position and {t) describes 
its displacement. 
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(39) 



V. N-BODY RESULTS WITH A MOVING 
STRAIGHT STRING IN EXPANDING 
SPACE-TIME 

We now turn to the cosmologically more relevant case 
of uniformly expanding space-time. Specifically, we con- 
sider an infinite, straight string moving through an ini- 
tially homogeneous distribution of particles expanding 
with the Hubble flow and follow the subsequent evolu- 
tion of those particles. 

We find that the energy momentum tensor T^"^ of the 
Nambu Goto string in a FLRW space-time is the same as 
in Minkowski space-time, except for a factor a~^. This 
factor cancels out when calculating T^^, which is used in 
the Einstein equations (|llll2p so that 



T, 



FLRW 



(37) 



It is straightforward to include the extra terms from the 
expansion that are present in the general equation (|16P '^ 
in the explicit equations of motion p3l34p . 

While the logarithmically divergent contributions at 
small and large distances in the metric perturbations 
(P7| and (pS)) did not enter the equations of motion in 
Minkowski space-time, we now have to deal with the term 
(d/a)Ei. We discuss in Appendix El the (unphysical) ori- 
gin of these divergences and how we regularise them, and 
we show that they do not infiuence the results. 



* a;* are comoving coordinates and the time derivative is w.r.t phys- 
ical time 



(at late times, a ^ ai, only the first term in the square 
brackets is relevant). We set e(j/o) = 1 for yo > and 
^ivo) = —1 for ?/o < as the velocity kick is always to- 
wards the string, Ui is the predicted velocity kick ((55)) 
and ti, tti are the time and scale factor at which the 
predicted velocity kick occurs (i.e. the moment when 
a^string = a;particie for a particular particle). From the par- 
ticle trajectory y[t) it easy to compute the turn-around 
radius, since at turn-around y = Q. One finds that for 
the particle that turns around at time t, 



(40) 



For testing purposes our simulation starts with a ho- 
mogeneous distribution of 32'^ particles on an uniform 
grid in a box size of L = 2Mpc//i. We start at redshift 
z = 10^. The string is initially at position (— 10L,L/2). 
We evolve the simulation with the string until it reaches 
position (10L,L/2) at which point we turn it off. There- 
after we evolve only the particles and since we use these 
simulations to test the implementation of our formalism, 
we use the expansion rate of a matter dominated universe 
throughout as the approximation above was derived un- 
der this assumption. 

Our results and comparison to the analytically- 
predicted trajectories and turnaround radii rt are shown 
in Figure [5j The numerically calculated particle trajec- 
tories initially follow closely the analytical predictions. 
Eventually the numerical results diverge from the ana- 
lytical ones due to the fact that in the simulation the 
velocity kick for particles that are further away from the 
axis of symmetry is smaller than the analytically pre- 
dicted velocity kick, which is derived assuming that the 
string that came from and went to infinity. 
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FIG. 3; Motion of a test particle in Minkowski space-time as the string passes by at a constant velocity v, from the adapted 
N-body code results. Initial conditions resemble the near by limit case, ie. the initial string particle separation is large in 
the X direction and small in the y direction. From the top to the bottom the panels show the particle position, velocity and 
acceleration. The left-hand panels show the x-component (along the direction of motion of the string) and the right-hand 
panels the y-component (perpendicular to the direction of motion of the string). 



VI. CONCLUSION AND OUTLOOK 

In this paper we set up the formalism needed to include 
topological defects and other sources of weak gravita- 
tional fields in N-body simulations. For this purpose, we 
have derived the equations of motion of massive particles 
(equation [15]) in a perturbed background. We find that 
both scalar and vector (gravito-magnetic) perturbations 
contribute significantly to leading order, {v/cf. Tensor 
perturbations on the other hand do not contribute to 
leading order, and their contribution to the acceleration 
is suppressed by one power of v/c. 

To test and illustrate the equations, we have applied 
the leading order result (equation [T5)) to the well known 
example of a straight string and we have recovered the 
expected velocity kick in the limit where the string passes 



very close to a particle (relative to the distance that the 
string travels). For the long, straight string we have also 
found that the vector contribution to the final particle 
motion is sub-dominant. 

Our results (eauations[T51and [T5)) can be used for calcu- 
lating the effect of any sources of weak gravitational fields 
on non rclativistic and on relativistic particles. The fo- 
cus of our continuing work is on doing large scale N-body 
simulations and calculating the effect of Abelian Higgs 
cosmic string networks (ssj on large scale structures. 
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Appendix A: Unphysical log(r^/ro) term 



An additional complication in the case of expanding 
background compared to the Minkowski space-time is 
evidenced by equation (jl6p . Since now d ^ 0, there 
is an additional force term proportional to E^, which 
is logarithmically-diverging with distance to the string. 
This unphysical divergence arises because of the idealised 
nature of the infinitely long, straight Nambu-Goto string: 
A string of finite thickness would regularise automatically 
the divergence for small radii, x"^ + y"^ ^ 0, while the di- 
vergence at large distances, x'^ + y'^ ^ oo is due to the 
assumed infinite length of the string. This is analogous to 
the logarithmic divergence exhibited by the electrostatic 
potential of an infinite line charge. The divergence at 
small separations is not an issue here because the physi- 
cal thickness of a realistic string is very small compared 
to the typical inter-particle separation in a cosmologi- 
cal N-body simulation. However, we would expect that 
in a cosmological setting the size of the causal horizon 
would provide an upper cutoff, above which strings are 
not formed. For this reason we have decided to regularise 
the logarithmic divergence in analogy with the potential 
of a line charge that is at a distance (P = x"^ + of a 
point charge and extends for a distance £ on either side, 
for which the logarithmic term becomes 



FIG. 5: Particle trajectories vs. scale factor. The solid lines 
are the simulated particle trajectories for several sample par- 
ticles at different initial distances from the passing string. 
The dashed lines show the corresponding trajectories based 
on the analytical predictions. The thick dash dotted line is 
the analytically-predicted turnaround radius. The analytic 
trajectories approach the axis of symmetry faster due to the 
fact that analytically all particles are treated as obeying the 
near by approximation. Their velocity kick is hence larger 
than the one calculated in the simulation. Furthermore, in 
the simulation we can see the second turnaround as the par- 
ticles fall back into the wake. 



Note also that for d ^ ^ the expression in brackets es- 
sentially is \og[cP / {2i)'^] which means that the integration 
constant in equations (P?]) and effectively corre- 
sponds to the string length 2£ in this limit. This term 
only affects the motion in the x direction, and hence it 
does not influence our results and comparison to the an- 
alytical predictions discussed above. 
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